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ABSTRACT 


A  number  of  tests  are  compared  for  testing  the  hypothesis 
of  a  constant  intensity  against  the  alternative  of  an  increasing 
intensity  function  in  a  nonhomogeneous  Poisson  process^-  The 
powers  of  the  tests  are  deteripined  by  Monte  Carlo  simulation 
against  alternatives  which  are  increasing  at  an  exponential 
rate,  a  power  rate  (Weibull  intensity],  and  a  logarithmic 
rate.  A  few  exact  powers  are  also  obtained. 

The  study  includes  the  well  known  Laplace  test  statistic 
which  is  known  to  be  appropriate  against  exponentially  increasing 
alternatives,  the  most  powerful  test  for  the  shape  ]Viramcter 
in  a  Poisson  process  with  Weibull  intensity,  the  likelihood 
ratio  test  against  arbitrary  NilJM’  al  ternat  i  ves  ,  two  non))a  ramet  r' i  c 
tests  for  trends  based  on  Kendall's  tau  and  Spearman's  rho, 
and  a  test  based  on  an  F-statistic. 
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1 .  INTRODUCTION 

Non-homogencous  Poisson  processes  (NHPPJ  ])rovide  models 
for  a  variety  of  physical  phenomena,  for  instance,  if  at 
each  failure,  a  system  is  repaired  to  its  condition  at  the 
time  of  failure  and  placed  in  service  again,  then  the 
failures  are  often  modelled  by  a  NIIPP  provided  the  rei)ai! 
times  can  be  neglected.  In  some  of  these  situations,  it  may 
be  reasonable  to  assume  that  the  intensity,  a(.),  is 
nondecreasing  and  so  tests  of  H^:  A(-)  is  constant  versus 
:  A(»)  is  increasing  are  of  interest.  This  would  indicate 
whether  the  simple  homogeneous  Poisson  process  (lIPP)  may  he 
adequate,  or  whether  a  more  general  NIIPP  model  is  required. 

Suppose  that  a  NHPP  is  observed  for  T*  units  of  time  with 

n  failures  and  failure  times  0  <  T,  I'  ^  ...  ^  T  1*. 

12  n 

(Of  course,  the  number  of  failures  is  a  random  variable,  which 
we  denote  by  N.)  We  compare  several  of  the  tests  which  aic 
available  in  the  literature  for  versus  in  this  time 
truncated  framework,  and  similar  conclusions  should  hold  for 
their  counterparts  based  on  data  truncated  on  the  number  of 
failures.  One  of  the  earliest  tests  is  attributed  to  haplace 
and  is  based  on  the  statistic  L  =  . T . /T* .  Under  IP, ,  the 

IP/T*  are  distributed  as  the  order  statistics  from  a  uniform 
distribution  on  (0,1),  so  L  behaves  as  the  sum  o f  un i fo rm 
random  variables  and  in  particular,  has  an  approximate 
normal  distribution  with  mean  n/2  and  variance  n/12.  Ilc'ntc-, 

H,.  is  rejected  if  L  >  n/2  +  z,  (n/12)^^^,  where  z,  is  the 


1-a  th  quantile  of  a  standard  normal  distribution.  Cox  (1955) 
showed  that  this  test  is  appropriate  for  testing  =  !)  versus 
3  ^  0  in  A(.t)  =  ae^^.  Bartholomew  (1  956)  gai^e  an  expression 
for  its  power  (see  also  Bates  (1955))  and  showed  that  it 
compares  favorably  with  the  one-sided  Kolmoj>orov-.Smi  rnov  test. 
Ascher  and  I'eingold  (  1978)  further  discuss  its  use  in  the 
study  of  repairable  systems. 

Another  family  of  intensity  functions,  wiiich  is  quite 
flexible,  is  A(t)  =  (3/6)  (t/0)^^  ^  for  3,0  -  0.  Because  this 
is  the  failure  rate  for  the  Weibull  distribution,  tlie 
corresponding  process  has  been  called  the  IVcibulI  I’oissois 
Process  (WPP).  Inferences  for  a  WPP  are  discussed  in  ('rou 
(1974  ,  1982),  Saw  (  1975),  Pinlcel  stc  in  (197b),  Lee  and  l,ee 
(1978),  Engelhard!  and  Bain  (1978)  and  Bain  and  laigclhardt 
(1980).  The  earlier  paper  by  Crow  gives  tests  for  e  with  e 
a  nuisance  parameter.  In  testing  3=1  versus  I-)  1  ,  wiiich  is 
equivalent  to  testing  versus  H^,  3=1  is  rejected  for 
small  values  of  Z  =  ^log  (T*/T| )  ,  which  has  a  chi-s()uared 

distribution  with  2n  degrees  of  freedom  under  and  is  lIMiTI 
in  this  WPP  setting. 

Boswell  (1966)  developed  the  likelihood  ratio  test  (I.RTI, 
conditional  on  n  failures,  for  11^^  versus  llj  and  an  ariiitrar\ 
NllPP.  The  maximum  likelihood  estimator  of  A(')  under  11,,^'  Il| 
is  shown  to  be  zero  on  fO,Tj)  and  constant  on 
k  =  1,2,. ..,n  with  =  T*  and 


A  computation  algorithm  is  also  given  there  (cf.  j).  1S()7). 

The  LRT  rejects  for  large  values  of 

W  =  ZfljJ^jlogCACTj.)  )  +  log(T*/M)  I  , 

2 

and  letting  x  denote  a  chi-squareJ  variable  with  k  degrees 

of  freedom 

lMW>w]  =  l'jJ^gHk,n)IMxhkM)  w| 

for  large  n,  where  the  P(k,n3  are  given  in  iable  A.S  o!‘ 

Barlow  et  al .  (1 972) . 


If  the  intensity  is  increasing  the  inter- fa i 1 ure  times, 

^  ^  l,2,.*.,n  (T^-O),  should  tend  to  decrease. 

Hence,  nonparamet ric  tests  for  trends  using  either  KendalJVs 
tau  or  Spearman’s  rho  could  be  considered  for  testing  11^^ 
versus  Hj.  (For  a  discussion  of  these  tests  as  tests  of  trend, 
see  Hollander  and  Wolfe  (1973,  p.  19(1).) 

Barlow  ct  al.  (1972,  p.  197)  observed  that  the  failure 
times  could  be  divided  into  two  parts  and  the  ratio 
(n“d)Tj/ (d  (T^-Tj) )  used  as  a  test  statistic  for  trend.  Oi' 
course,  an  increasing  intensity  should  correspond  to  a  larger' 
value  of  this  ratio.  If  the  intensity  is  constant ,  then 
conditional  on  n  failures,  T,  can  be  expressed  as  ,  V  .  / |  V  . 
where  the  Yj  arc  independent  exponential  variables  with  a 
common  mean.  Hence,  is  rejected  if  tiris  ratio,  which  we 
denote  by  F,  exceeds  the  1-m  th  quantile  of  the  I*  distribution 
with  2d  and  2(n-d)  degrees  of  freedom.  We  consiilered  t  lu'  tc'st 
with  d  =  [n/2]  and  the  one  with  Tj  ^  T*/2  •  However,  the 
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latter  test  is  only  applicable  if  both  of  the  intervals 
(0,T*/2]  and  (T*/2,T*]  contain  failures  and  the  former  test  is 
applicable  whenever  n  ^  2.  For  small  truncation  times  this 
difference  can  be  appreciable  and  so  in  that  which  follows, 
we  only  consider  d  =  [n/2]. 

Section  2  contains  the  results  of  a  study  which  compares 
the  powers  of  the  tests  described  above  and  in  Section  S, 
our  recommendations  based  on  this  study  are  given. 


POWIiR  COMPARISONS 


Some  of  the  tests  discussed  in  the  Inst  section  are 
clearly  appropriate  for  certain  families  of  intensities, 
for  instance  the  test  based  on  Z  is  IJMPl)  in  the  WPP  setting, 
and  it  would  be  of  interest  to  study  their  j^ower  functions 
for  a  broad  range  of  intensities.  On  the  other  hand,  soiik'  oi 
the  tests,  such  as  W,  are  designed  to  be  more  oitinihu>  and 
information  concerning  the  loss  of  power  incur  ied  uhen  us  ins, 
them  instead  of  a  test  which  is  optimal  in  a  particular 
situation  would  be  useful. 

Because  of  the  complex  nature  of  some  ol*  the  powcM*  f'unction 
involved  a  Monte  Carlo  study  was  conducted.  I'lie  tests  based 
on  L,Z,W,F  Kendall's  tau  (which  we  denote  by  K)  and  Spearman's 
rho  (which  we  denote  by  S)  are  compared.  fhe  following,  lesull 
is  the  key  to  generating  a  NUPP  with  intensity  A(-).  If 

rt 

A(t)  =  A(s)ds  is  strictly  increasing  and  (Su  1  denotes  the 
failure  times  for  a  HPI'  with  A  =  l,  then  {A 

failure  times  for  NHPP  with  intensity  •  )  .  (‘(t)  is  the 

expected  number  of  failures  in  ( 0 ,  t ) . )  Intensity  functions 
which  are  exponential,  of  tiie  Weibull  tyjK',  logarithmic  and 
those  that  increase  but  have  a  horizontal  asymptote  are 
considered.  A  particular  value  of  T*  may  be  appropriate 
for  one  intensity  function  but  not  for  another,  that  is  it  ma> 
give  rise  to  a  reasonable  amount  of  data  I'or  one  intensit\ 
function  and  no  data  (or  a  very  large  sample  size)  for  another. 
So  the  truncation  time  will  be  allowed  to  vary  with  M-)  and 


in  particular  T*  is  chosen  so  that  the  expected  sample  si::e, 
h(N) ,  is  10,  20  or  40. 


If  A[*)  is  an  intensity  function  with  correspond  i  ni^  mean 
function,  A(')  and  X^(t)  =  A(t/0)/(3,  then  A^^(t)  =  A(i't).  Hence, 
the  failure  times  for  the  process  with  intensity  can 

be  expressed  as  A  ^(S^)/0  where  the  S.  are  the  failure  times 
for  a  Poisson  process  with  A=l,  All  of  the  tests  consideied 


are  scale  invariant  and  so  we  need  onlv  consider  (;=J.  I  j) 


particular ,  we  consider  processes  with  intensities  of  the 
form  A(t)  =  iSt^  ^  with  B=l,2,4,  A(t)  =  ae^  with  .A, 1,2 
A(t)  =  (x^  V  (ot- 1 )  ! )  /  ^  j  /  j  !  )  ,  which  is  the  failure  rate 

of  a  gamma  distribution,  with  nt=2,4  and  A(tj  =  log(t  +  l). 

The  gamma  intensity  was  chosen  because  it  increases  slowly  and 
in  fact  approaches  a  constant  as  t  flowcver,  for  such 

alternatives  the  data  beyond  some  point  in  time  resembles 
that  for  a  constant  intensity.  The  tests  considered  do  not 
have  enough  power  in  such  cases  to  make  interesting  comparisons. 
If  one  wishes  to  discriminate  against  this  kind  of  an  intensity 
several  replications  of  the  process  should  be  observed  for  an 
initial  time  interval  rather  than  observing  one  process . 

For  this  reason  the  study  of  the  GPP  was  not  completed  and  no 
further  results  concerning  it  are  given  here . 

Table  1  contains  Monte  Carlo  estimates  of  the  powers  for 
the  tests  and  alternatives  described  above.  These  estimates 


are  based  on  5000  iterations.  Tests  with  a  nominal  significance 
level  of  a  =  .05  are  considered,  however  the  LRT  is  known  to 
have  larger  true  signi f icance  le\ el .  For  instance,  wi  r h 


with  HtN)  =  10  and  a  target  level  of  .05  its  esti/nated  significance 
level  is  .0()4.  bxper  iment  i  ng  we  found  t  hat  a  noiiii  na  1  level  ol 
N'ields  an  estimated  level  of  a[)prox  imatcl  y  .05  for  the  idfi  and  the 
range  of  sample  sizes  being  considered.  Clearly,  this  is  a 
disadvantage  of  the  LR'f  in  this  situation.  As  can  be  seen  from 
Table  1  the  significance  levels  for  the  other  tests  arc  very 
close  to  .05.  Of  course,  the  Z  and  f  tests  are  exact. 

To  give  some  idea  of  the  accuracy  in  these  est imatcul 
powers,  we  consider  the  power  of  Z  for  f  >  1.  Since 
A(t)  =  t^,  the  failure  times  have  the  same  joint  distributions 
as  with  the  S.  as  before,  and  T*  =  •  Conditional 

on  n  failures,  the  distribution  of  Z  is  the  same  as  Z/i  under 
\(t)  1.  So  the  power  of  Z  at  >:  is 

(->,“(  2n  )  |/(n)  )  )  , 

2 

where  l^he  a-th  quantile  of  a  chi-squared  distribution 

with  2n  degrees  of  freedom.  The  first  row  in  Table  1  gives  the 
cAact  powers  computed  from  this  formula.  I’hc  largest  discrepancy 
between  these  exact  and  estimated  powers  is  .004. 

Similar  comparisons  can  be  made  for  the  L  test  since  it  is 
based  on  approximate  normal  distributions.  Yor  a  NlIPP,  given 
that  there  are  n  failures  in  the  interval  (0,T^],  the  failure 
times  are  distributed  as  the  order  statistics  of  a  raiuloni 
sample  of  size  n  from  the  density  f(t)  =  A(t)/A(T*),  0  •  t  ■  T* 
and  f(t)  =  0  otherwise.  Since  I.  can  be  expressed  as  the  sum 
of  the  unordered  observations  from  this  density  divided  by  T*, 
it  has  conditional  mean 

fT^ 

U  -  i;(L|N  =  n)  =  n  t  A  (t)dt/ tT*A  (T^l) 

n  J,) 

I 


and  conditional  variance 

o“  =  V(l.lN  =  n)  =  n{  t  ^  (  t )  d  t  /  ( T*  (T*  )  )  -  u  “  I  . 
n  j  n 

Hence,  the  power  of  L  is  approximately 

(A(T*33"{l-n(n/2-p^^  +  =  j_  jn/12)^/"j/-  :  /  (  ii  I  (  1 -e'  ' 

2 

However,  in  the  WPP  case  the  tormulas  for  ;i  and  o  siini')lir\'  to 
’  n  n  ^  ' 

=  nB/(^'''^l)  and  =  n  [  B/ (  B+ 2) -(.:/(  1  J  J  ^  ]  . 

The  fourth  row  in  Table  1,  which  is  labeled  L  ^appi'ox.),  yives 
the  approximate  j^owers  comj'JUted  using  the  formulas  and  the  fil'tli 
row  (labeled  L(est.))  gives  the  Monte  Carlo  estimates  oi' 
these  powers.  The  largest  discrciKincy  !')etween  the  a[)ju-ox  i  iiia  t  c' 
and  Monte  Carlo  values  occurs  for  b(N)  =  10  and  B=2  and  the 
difference  in  this  case  seems  to  be  due  primarily  to  the  normal 
approximation.  One  could  obtain  a  further  check  on  accurac)'. 
Bartholomew  (1955)  gives  the  power  of  L  for  exponential 
intensities,  conditional  on  n  failures,  and  this  could  be  used 
with  the  law  of  total  probability  to  obtain  a  formula  for  the 
power  of  L  for  such  alternatives.  Because  of  the  comj’^lex  nature 
of  the  conditional  power  this  has  not  been  carried  out. 

To  aid  in  the  comparison  of  these  tests,  relative  elficiencie 
are  estimated.  For  each  case  considered,  for  instance  the 
WPP  with  3  =  2  and  F(N)  =  10  is  one  case,  the  relative  efficiencx’ 
of  a  test  is  its  power  divided  by  the  largest  of  the  six 
powers  for  that  case.  Table  2  contains  the  estimated  relative 
efficiencies,  that  is  the  ratio  of  the  estimated  powers.  The  Z 
test  was  the  most  powerful  for  the  loga r 1 1 hmi c  and  We i bu 1 1 


( 1 II  the  ca  se  IVPP  , 


aiuJ  1.  (  X  ) 


20  tlie  Tafl 


a  1 ternat  i  ves  . 


that  the  estimated  relative  crficieney  of  Z  is  less  tfiaii  one 
is  due  to  Monte  C’arlo  error.  Recall  Z  is  IIMIMI  in  tlii>> 
setting;  ‘And  the  true  power  is  .oopo.)  The  L  test  uas  most 
powerful  for  the  e\{)onential  intensities. 

It  should  be  noted  that  the  mininium  of  the  relative 
efficiencies  of  Z  in  our  study  is  .014,  (We  denote  niiniiiium 
relative  efficiency  by  MRh.J  'I'his  is  the  la  ripest  MRl.  of 
the  six  tests  studied  and  so  we  recommend  the  Z  test  if  oi)c‘ 
wishes  to  discriminate  against  intensities  in  the  I'ange 
considered.  Ihe  MRI-  for  the  hR'i’lW)  is  .810,  for  h  it  is  .8Zo, 
for  Spearman’s  I'ho  (S)  it  is  .530  and  for  Kendall’s  tau  (K) 
it  is  .521  and  for*  t(ie  f  test  it  is  .  522  .  If  one  is  not 
concerned  with  the  possibility  of  slowly  increasing  intensities, 
such  as  the  logarithmic  intensity,  then  the  h  test  could  be 
used  which  has  an  estimated  MRI:  over  that  range  of  .05.  llo\\ever, 
the  increase  in  MRh  is  not  large  and  it  seems  to  be  best  to  use 
Z  and  protect  against  the  possibility  of  slowly  increasing 
intensities.  Vor  the  situations  considered  here  the  non]ni  ranie  t  r  i  c 
tests  can  not  be  recommended,  but  if  the  distribution  assuiiijUion 
on  inter- fai  lure  times  is  in  question,  one  might  wish  to 
consider  one  of  them.  (There  is  very  little  difi'erence  in 
their  MRl:s.)  The  1*  test  performs  better  for  large  h(N)  tlian 
for  smaller  H(N),  but  its  performance  for  slowl)*  increasing 
intensities  is  not  strong.  The  LRT  was  designed  to  discriminate 
against  all  nondec rcas i ng  intensities  which  are  not  constant. 

If  one  were  concerned  about  very  nonregular  intensities  W  could  be 
considered,  but  such  intensities  were  not  considered  here. 


CONCl.USlONS 


3 . 

In  testing  for  an  increasing  intensity  in  a  Poisson 
process,  the  1  test  performs  quite  well  for  the  range  of 
alternatives  studied  here,  that  is  for  logarithmic,  Weilnill 
and  exponential  intensities.  In  fact,  its  efficiency  relative 
to  the  other  five  tests  considered  is  at  least  POu  and  for 
the  logarithmic  and  Wcibull  intensities  it  is  the  most 
powerful  of  the  six.  I’or  alternatives  of  the  type  studied 
here  the  Z  test  is  recommended. 
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TABLli  1.  Estimated  powers  for  testing  versus  Hj 
Weibull  Intensity:  A(t)  =  Bt^  ^ 


I.O 

B 

2.0 

4 . 0 

Test\. 

10 

2  0 

40 

10 

20 

4  0 

10 

2  0 

4  0 

Z (exact) 

.  0500 

.  0500 

.  0500 

.6209 

.  8989 

.  9957 

.  978  5 

.  9999 

I  .  OOOD 

Z(est.) 

.  0S2 

.  053 

.052 

.  624 

.902 

.  99  5 

.  7  () 

.  9  96 

1  .  ()()() 

iv^ 

.  052 

.048 

.050 

.  511 

.797 

.  978 

.953 

.  9  9  9 

1  .  000 

L(approx.) 

.  0500 

.  0500 

.  0500 

.  5721 

.  8558 

.  98  93 

.  9787 

.  9  9  9  9 

1  .  0000 

L (est . ) 

.  054 

.050 

.  047 

.  593 

.855 

.  98  9 

.  97  3 

1  .  000 

1  .  0  00 

S 

.  049 

.050 

.049 

.331 

.534 

.783 

.  542 

.774 

.958 

K 

.  046 

.  052 

.052 

.325 

.  539 

.  786 

.  541 

.  78  2 

.  957 

1- 

.  052 

.  046 

.046 

.  381 

.614 

.  905 

.852 

.  9  90 

1  .  00(1 

Exponential 

Intensity:  A{t)  =  uc 

(1 

t 

0.5 

1.0 

2.0 

\E(N) 

Test"'  ' 

10 

20 

40 

10 

2  0 

40 

10 

20 

4  0 

Z 

.  762 

.981 

.997 

.608 

.  940 

.  996 

.4  36 

.844 

.  99  4 

.729 

.984 

1.000 

.  582 

.939 

1 . 000 

.400 

.818 

.998 

L 

.805 

.993 

1.000 

.656 

.972 

1.000 

.477 

.8  92 

1 . 000 

S 

.484 

.  832 

.  986 

.395 

.778 

.  981 

.2  90 

.670 

.  96:- 

K 

.476 

.833 

.987 

.383 

.777 

.981 

.274 

.  668 

.966 

F 

.620 

.956 

1.000 

.481 

.892 

1.000 

.343 

.  74  6 

.  99  5 

Logarithmic  Intensity 

X(t)  =  log(t+l) 

•--.^.__l’es  t 

1 

E(N)^ 

Z 

W* 

1, 

S 

K 

I' 

10 

.  562 

.311 

.321 

.  205 

.202 

.206 

20 

.529 

.444 

.439 

.2  94 

.  298 

.2  76 

40 

.770 

.650 

.643 

.426 

.429 

.422 

1.  The  nominal 

level 

for  the 

W  test  was  .04 

« 

TABLli 


1 


Hstimated  relative  e f f i c i ene i es 


Log, 

arithmic 

Intensity 

Wc i hul  1  Intensity: 

'  (  t  )  -  t 

,  f-  -  1 

A(t)=log(t+l} 

2.0 

4 . 0 

H(N) 

rciTt^ 

10 

2  0 

40 

10 

20 

4  0 

10 

2  0 

40 

'7 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1.000 

1 . 000 

1 . 000 

.  9  96 

1  .  000 

.859 

.  859 

.844 

.819 

.884 

.  98  3 

.  976 

.  9  9  9 

1 . 000 

L 

.888 

.8  50 

.835 

.  950 

,  948 

.  994 

.  997 

1  .  000 

1.00  0 

S 

.  56b 

.5  56 

.  553 

.  530 

.  592 

.  787 

.555 

.  :’74 

.  9  58 

K 

.  558 

.  565 

.  557 

.  521 

.  598 

.  790 

.  554 

.  782 

.  9  5  7 

{■ 

.  569 

.  52  2 

.  548 

.611 

.  681 

.  9 1  0 

.  875 

.  990 

I  .  000 

Lxponent ial 

Intensity:  A ( t }  =  ae 

t 

0.5 

a 

1 . 0 

2.0 

h(N) 

Test 

10 

20 

40 

10 

20 

4  0 

10 

2  0 

4  0 

Z 

.947 

.988 

.997 

.927 

.  967 

.  996 

.914 

.  946 

.  994 

.906 

.  991 

1.000 

.887 

.964 

1.000 

.8  39 

.917 

.  998 

L 

1.000 

1.000 

1.000 

1.000 

1.000 

1 . 000 

1.000 

1  .  000 

1 . 000 

S 

.601 

.838 

.  986 

.602 

.800 

.  981 

.608 

.751 

.  Ob” 

K 

.  591 

.  839 

.987 

.  584 

.  799 

.981 

.  5  74 

.  749 

.  96  6 

I- 

.770 

.  963 

1.000 

.753 

.918 

1 . 00  0 

.719 

.8  56 

.  99  5 

^The  nominal  level  for  the  W  test  was  .04. 


